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Contractive states for a free quantum particle were introduced by Yuen (Yuen H P 1983 
Phys. Rev. Lett. 51, 719) in an attempt to evade the standard quantum limit for repeated 
position measurements. We show how appropriate families of two- and three component 
\ "Schrodinger cat states" are able to support non-trivial correlations between the position 

£Nj ■ and momentum observables leading to contractive behavior. The existence of contractive 

^ ' Schrodinger cat states is suggestive of potential novel roles of non-classical states for precision 

measurement schemes. 
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I. INTRODUCTION 

m 
o , 

fT) \ Shortly after the seminal paper by Schrodinger in 1935 lj, "Schrodinger cat" became a pictorial 

way to refer to a prototypical family of genuinely non-classical states i.e., quantum states without 
a classical counterpart. While the introduction of such states found its original motivation in the 
analysis of the celebrated cat Gedankenexperiment in quantum measurement theory 0, 0], the 
quantum-mechanical weirdness of Schrodinger cat states has since then been widely appreciated in 
various contexts. The latter range from quantum non-locality ^3], to non-classical states of light 
and matter ;4], and to the emergence of classicality from the quantum world [H, E|. in particular, 
cat states arising from the superposition of macroscopically distinguishable states of a quantum 
degree of freedom play a prominent role in the debate on macroscopic quantum mechanics 0]. 
More recently, the rapid development in the field of quantum information science 0] has prompted 
the consideration of various classes of non-classical states from the point of view of the quantum 
resources that they may involve • EH- For instance, cat-like superpositions of bosonic coherent 
states have been shown to provide robust encodings against amplitude damping ^3], to carry the 
potential for improving the sensitivity of weak-force detection [131 ] . as well as to exhibit generalized 
entanglement properties 

While applications of non-classical states within quantum information processing are still emerg- 
ing, their association with the field of high precision measurements has a long history. A leading 
motivation stems from the need for understanding the ultimate resolution limitations for repeated 
position measurements in experimental and observational gravitation Without a careful prepa- 
ration of the initial state, quantum mechanics sets a limit based on the assignment of equal un- 
certainties to position and momentum, resulting in the so-called Standard Quantum Limit (SQL) 
[ lfi \vfy . Attempts to beat the SQL have stimulated the development of a theory for quantum 
non-demolition measurements LUJ], and several experimental proposals aimed at improving the 
sensitivity of resonant detectors of gravitational radiation j^. Efforts have primarily focused on 
a harmonic oscillator ^J], resulting in the possibility of overcoming the SQL using non-classical 
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squeezed states. The case of a free test mass did not receive much attention until Yuen proposed a 
novel class of so-called contractive states analogous to the bosonic two-photon coherent states |2lj ] 
and capable, in principle, of narrowing the position variance below the SQL bound. While Yuen's 
proposal initiated a controversial debate 23, 3, 3, iH Ef3, 22], 2^1 , the demand for schemes able 



to beat the SQL for free masses has meanwhile increased due to the amazing improvements in the 
sensitivity of gravitational wave interferometric detectors |2j|. In principle, even a modest gain in 
the sensitivity obtainable by exploiting non-classical states would imply a relatively large increase 
in the fiducial volume for gravitational wave detection [iO, 0, Unfortunately, contractive 

states as proposed by Yuen are difficult to produce in practice; the only concrete proposal that has 
appeared so far in the literature was regarding their generation for atomic states 33]. 

In this paper, we revisit Schrodinger cat states from the motivating perspective of the SQL of 
a free mass, and show how in addition to their already known properties they are also capable of 
manifesting contractive features. The content is organized as follows. In Section II we summarize 
the basic notions about the contractivity property as introduced by Yuen. In Section III we 
investigate the behavior of the position variance for a paradigmatic class of Schrodinger cat states, 
evidencing its dependence upon various parameters, and then outlining possible directions for 
generalization of this result. Section IV is devoted to discussing some implications of the existence 
of contractive Schrodinger cat states, and preliminarily assessing the prospects for demonstrating 
contractive cat behavior. Additional remarks are presented in the conclusions. 



II. CONTRACTIVITY AND YUEN STATES 

Consider a quantum particle of mass m freely evolving in one dimension. If x,p are the position 
and momentum observables, the equation of motion for x in the Heisenberg picture is given by 

x(t)=x(0) + ^t, t>0. (1) 
m 

Let Ax 2 and Ap 2 denote the variances of the corresponding operators i.e., Ax 2 = {Ax 2 ), Ap 2 = 
(Ap 2 ), with Ax = x — (x), Ap = p — (p), and (x), (p) being the average position and momentum 
values, respectively. According to Eq. ([T)l. the position variance after a time t can be related to 
the initial uncertainties Ax 2 (0) and Ap 2 (0) by 

Ax 2 (t) = Ax 2 (0) + 1 2 > Ax 2 (0) , (2) 

expressing the fact that the initial position variance increases under free evolution. The SQL 
for position measurements coincides with the minimum value attainable by Ax 2 (t) in (J2J) when 
t is regarded as the time interval between two successive, identical position measurements j^. 
The argument is based on a heuristic application of the Heisenberg uncertainty principle, whereby 
Ax(0) is interpreted as the uncertainty due to the finite resolution of the first measurement, and 
Ap(0) represents the momentum disturbance caused by the measurement (see, however, [2^] for a 
substantially more refined treatment). Because Ax(0)Ap(0) > h/2, Eq. can be rewritten as 

h 2 t 2 

Ax 2 (t) > Ax 2 (0) + - OA . , (3) 
w - v; 4m 2 Ax 2 (0) v ' 

which, upon minimization with respect to Ax 2 (0), implies 

Ax 2 (t) >^ = Ax 2 QL (t) . (4) 
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Accordingly, the outcome of the second position measurement cannot be predicted with uncertainty 
smaller than the SQL value, Ax(t)sQL = (ht/m) 1 ^ 2 . 

Yuen pointed out a major flaw in the above derivation, as by adding in quadrature the initial 
position variance and the one propagated in time from the initial momentum variance as in Eq. 
(J2J) it is implicitly assumed that position and momentum are initially uncorrelated. While this is 
correct, for instance, in the common case where the mass starts in a minimum-uncertainty state, 
it is necessary in general to relax this assumption, by replacing (J2J) with the full expression: 

Ax 2 (t) = (m\m 2 \m) - (m\Ht)\m) 2 

= Ax 2 (0) + 

m m z 

where { , } denotes the anti-commutator operation. Thus, the position variance is directly sensitive 
to the initial correlation coefficient A5L(0), with 

A 2 p (0) = ({Ax(0),Ap(0)}} = ({x,p})-2(x)(p) = -ih + 2(xp) - 2(x)(p) , (6) 

upon explicitly using the canonical commutation rules for x,p [3J|. Yuen also proposed an explicit 
class of quantum states, the so-called twisted coherent states, where a negative correlation term 
A 2 p (0) < is realized (21j. As a consequence of the parabolic dependence of Ax 2 (t) upon t implied 
by ©, the initial position variance then narrows, attaining a minimum at an optimal time i > 0. In 
addition, the variance remains below its initial value Ax 2 (0) for the finite time interval < t < 2i, 
opening the possibility to evade the SQL for appropriate values of the parameters. 

It is interesting to note that the essence of Yuen's proposal relies on allowing complex Gaussian 
states. In one dimension, for instance, consider a wave function of the form 

ip{x) = k exp (— ax 2 ) , (7) 

with a € C. By writing a = a% + ict2, ai,ct2 £ ^> the normalization condition leads to k = 
(2ai/7r) 1//4 (up to an irrelevant phase factor). Because for such a state {x) = 0, {p) = 0, and 
(xp) = iha/2ai, the corresponding correlation coefficient becomes 

A 2 xp (0) = h^^ = -2H, (8) 
1 ct\ 

where the parameter £ = a^l2a\ can be, in principle, arbitrarily large. Note that the value 
£ = corresponds to the standard Gaussian state, while £ > implies a negative correlation 
term responsible for the desired non-monotonic behavior of Ax 2 (t). Following Yuen's notation, the 
prototype wave function Q can be generalized to the following family of states: 

^ (x)= ( vr%-H 2 ) ' eXP {-^J^ iX - X ° )2 + lE i iX - X0) } ' (9) 

where now {x)y = xo, (p)y = Po are real numbers, and the parameters /i,^, and to are related to 
the initial position and momentum variances and correlation as follows: 

^■ {0) = ^-"^ (10) 

. 2 , n , hmuj 1 + 4£ 2 
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A^,y(0) = -2fi£ . (12) 

Note that a% = mui/2h\fi — v\ 2 (compare Eq. Q), and the position variance evolves explicitly in 
time as: 



h 2 2h i t+ fiw l + 4£ 2 

2muj 1 m 2m |/i — ^| 2 



Aa^r(t) = - l/u-^l i + - — -j %r. (13) 



Motivated by the comparison with the SQL, a natural strategy is to optimize the relative 
variance obtained upon normalizing Ax Y (t) by the spreading expected at time t for a state evolving 
at the SQL as given in (j3J) i.e., to minimize the quantity 

. - N Axl(t) Axtit) 

By differentiating Ay(i), we get a minimum relative variance at a time 

# |/x — z^| 2 , . 

V = u,(l + 4^2)1/2 • ^ 

This criterion coincides, for large values of £, with the direct optimization of the absolute variance 
(fT3l) originally adopted in [2l| : 

thanks to the fact that the ratio tyjiy = (1 + 4£ 2 ) 1/,2 /2£ quickly approaches unity with increasing 
£. The relative criterion (|14jl. however, has the advantage of allowing for a continuous comparison 
with the SQL value at a generic instant of time t. From equations ()13|) and (|15jl. the minimum 
relative variance is therefore 

A y (t*) = (l + 4e 2 ) 1/2 -2e, (17) 

which quantifies the advantage in the use of a Yuen contractive state with respect to one with 
zero (or positive) correlations. As a function of £, the ratio in (|17|) is monotonically decreasing 
in the contractivity region £ > 0, a unit value corresponding to a classical state with £ = 0. 
While arbitrarily small values can be attained asymptotically, a significant gain already appears 
for relatively small £ parameters; for instance Ay(ty) ~ 0.41 for £ = 1/2. 

As mentioned, Yuen's proposal sparked a lively debate on the meaning and conceptual consis- 



tency of the SQL for a free mass [22J, [23J, |2J, |25|, |26(], and the discussion still continues |27|, |2£ 
While various issues have been substantially elucidated, from an operational perspective limited 
effort has been devoted to address possible ways for preparing Yuen states in the laboratory. To the 
best of our knowledge, the only steps in this direction have been taken by Storey and co-workers 
[ssT ]. who proposed a position measurement scheme for atoms entering a strongly detuned standing 
light wave in a cavity, and by Vitali et al. in the context of stochastic cooling of macroscopic 
mirrors. Ultimately, one obstacle encountered in designing realizable preparation schemes for Yuen 
states is related to the difficulty of dealing with a complex a. Thus, although very interesting in 
principle, this class of states seems not very palatable to experimentalists, motivating the search 
for alternative families of contractive states. 
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III. CONTRACTIVE FEATURES OF SCHRODINGER CAT STATES 

Following Yuen's original suggestion, position-contractive states for a free mass are denned by 
the property 



= Aj„(0) < . (18) 

t=0 



Because of the equation of motion @, Ax 2 (t) increases monotonically with time whenever A^ p (0) > 
0. Equivalently, in terms of the relative variance A(i), one may verify that if A^ p (0) > 0, the 
inequality A(t) > 1 for all t is also enforced as a consequence of the Heisenberg uncertainty 
principle, and vice versa. In particular, the condition l|18|) is never satisfied by a preparation in 
a classical Gaussian state. However, Yuen's class is not the only family of non-classical Gaussian 
states one may a priori consider. Another notable class capable, as we shall describe now, of 
fulfilling (fTH|) is provided by macroscopically distinguishable Gaussian states i.e., Schrodinger cat 
states in the position variable. 



A. Two-component cat states 

Let us consider, to begin with, a two-component Schrodinger cat state, namely a coherent 
superposition of two one-dimensional Gaussian wave packets with inverse variance parameter a € 
M + , localized around two positions +xq and — xq, xq > 0, and having relative amplitudes k + , k- € 
C: 

^S2(x) = N 2 {k + exp[-a(x - x ) 2 ] + k- exp[-a(x + x ) 2 ]} . (19) 

Here, the normalization constant is given by 

\N 2 \ 2 = (2a/7r) 1 / 2 [\k + \ 2 + \k_\ 2 + exp(-2axl)(klk^ + k*_k + )y 1 , (20) 

and we have assumed that the two Gaussians have zero initial momentum. Note that, with this 
notation, each macroscopically distinguishable component in |V'52(^)| 2 corresponds to a Gaussian 
probability distribution with variance l/4a. After straightforward calculations we obtain the initial 
average values of position, momentum, and their product as follows: 

| k.^ j \k — | /c\t\ 
{X)S2 = \k+\ 2 + |fc_| 2 + exp(-2a4)(klk- + k*_k+f° ' { } 

(P)s2 = 4 ^ a | fc+ | 2 _ jfr" | 2 exp(-2axo)(x) 52 , (22) 

(xp) S2 = ih/2 . (23) 
Thus, the correlation coefficient for a mass initially prepared in the state (|19|) is: 

lm(k+k* ) 

K P ,S2(°) = -2<x)s2(j5)s2 = ~8to |fc+|2 V + |fc - ; |2 exp(-2axg)(x)| 2 . (24) 
According to (|24|). two contractivity regions fulfilling (|18jl are then possible in principle: 
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This illustrates a different route for generating a non-zero correlation coefficient compared to 
Yuen states. In the latter case, regardless of whether (x)y{P)y = or not, a non-zero correlation 
^xpY ar i ses due to the asymmetric contribution from (xp)y and (px)y, enforcing ({x,p})y ^ 0. 
For Schrodinger cat states of the form (|19|). it is instead the initial anti-correlation between x and 
p, {{x,p})s2 = 0, that induces contractivity for appropriate k + ,k- and (x)s2 7^ 0. From a less 
formal viewpoint, in the case of the Yuen states the presence of an imaginary part in the quadratic 
coefficient a effectively correlates the position and momentum observables. For cat states as in (|19[1. 
the correlation is obtained through the delocalized spatial structure, which can enforce non-zero 
average values for both x and p as in (|2*T]) - (|2*2*|) . 

In order to evaluate the position variance at time t based on ©, it is necessary to further 
calculate the second moments Ax| 2 (0) and Ap| 2 (0) ( see a ^ so the appendix). For convenience, we 
set a = 1/2 Aq henceforth, and introduce an adimensional time parameter as 

It is worth noting that while the values of xq and Ao separately affect the expectations (x)g2, 
(p)s2, the correlation coefficient (and hence contractivity) depends only upon their ratio, which 
we denote by 5 = xq/Aq. We can also assume, without loss of generality, that k+ € R. By letting 
k + = c+, k- = C-e~ ld , with c±, 6 € M and 9 determining the initial relative phase between the two 
Gaussian components, we finally reach the rewritten forms: 

„ 9 / s An 2c+c_ + (ci + c 2 )e~ 5 ' 2 cos 6* 9 , . 

Ax| 2 (0) = -f + 2c + c_^ D x ° ' ( 27 ) 



ApI 2 (0) = ^ 



A^ 2c + c-e~ 52 + (c\ + c 2 )cos6» _ S 2 2 
— 2c + c_ — e x 



(28) 



In the above equations, D 
tion variance is then: 



cL + c 2 _ + 2c+c_e &2 cos 



. The complete expression for the posi- 



AX S2(V) 



-7T- + 2C+C_ 



A 



2c + C- 



2c + c_ + (c\ + c 2 _)e~' 5 cos 
2c+c_e~ 52 + (c 2 ^ + c 2 .) cos# _ g2 



D 



e Xr 



4c+c 



v 2 - 



d - c 2 



<5 2 „ 2 



■ sm t/e x. 



(29) 



Thus, the variance may be written as Arc| 2 (ry) = A + Br] + C?7 2 , where the coefficients A,B, C 
of the relevant powers of r] can be inferred from the quantities in square brackets in Eq. (j29|) . 
By construction, A reproduces the initial position variance, while B and C are responsible for 
the spreading from the initial value. The behavior for a single initial Gaussian state may be 
recovered by simply setting either c+ or c_ to zero. Note that, regardless of possible contractive 
features, the rate of spreading of Ax? ?2 (?7) may be slower than in the Gaussian case due to the 
interference term contained in C. By normalizing Ax^if]) to the reference value dictated by the 
SQL, Ax 2 S Q L (rj) = AqT], we obtain: 



(30) 
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where A2 = A/Aq etc. By minimizing As2('n) with respect to 77, an optimal time equal to 77* = 
(A2/C2) 1 / 2 is found, which is independent of the linear coefficient B2. The corresponding minimum 
value attained by A$2 is 



A S2 (r/*) =B 2 + 2(A 2 C 2 ) 



1/2 



(31) 



which can be further optimized with respect to the remaining parameters. 

To analyse in detail the behavior of A$2, it is convenient to additionally introduce the ratio 
k = c+/c—, and rewrite the relative variance (|29j) as: 



As2(»7>K,0,£) 



2 2K + (1 + At 2 )e 



- + 2k5' - 

2 (1 + k? + 2ne- 52 



-S 2 



cos ( 



+ 



2k6' 



2ne- 



+ (1 + k 2 



COS ( 



1 COS I 



1 

- + 

V 



[1-K 2 )e 



-J 2 



sun 



(1 + k 2 + 2«e 



-<5 2 



cos ( 



(1 + K 2 + 2kc s cos 



(32) 



where the dependence of A52 upon the relevant parameters has been made explicit, and the ex- 
pressions in the square brackets now correspond to A2,B2,C2 of (|3Uj) . According to (|32j) . A52 is a 
complicated function of both the time 77 and the various parameters k, 5, 9 which characterize the 
initial state. In terms of the new parametrization, the contractivity regions given in (|25[) become: 



(I) 



K > 1 , 

sin 9 > 



or 



(II) 



K < 1 , 

sin 9 < 



(33) 



Note that the dependence upon 9 is periodic, thus we can limit the analysis to the range 
0° < 9 < 360°. In addition, A52 satisfies the following invariance property: 



A S2 (r], k, 9, 8) = A S2 (r 1 , k' 1 , 360° - 9, 5) 



(34) 



Thus, given the behavior of A^2 in a single contractivity region, say (I), the behavior in the remain- 



ing region is the same upon transforming k 



and 9 1— > 360° — 9. Some qualitative insights 



into the behavior of the function 1)32(1 can be obtained by inspection of the various terms. The 
contractive term, B2, assumes zero values at 9 = 0° and 180°, namely, when the two distinguishable 
components of the cat add in phase or anti-phase, as well as when k = 1, which corresponds to an 
equally weighted superposition state. As a function of 9 and 5, |B2 [ is large for 9 ~ 90° and 5 of 
the order of unity. 

By numerical analysis of Eq. (|32|). the figure of merit for contractivity A52 is found to attain 



its minimum 



A^ n (r]*, K *,9*,S*) - 0.757, 



(35) 



for optimal parameter values rj* ~ 1.105, k* ~ 2.26, 9* ~ 127°, and 5* ~ 0.49. The dependence 
of A52 upon the effective time 77 and k is displayed in Fig. Q for 9 = 9* and 5 = 5*. The 
section in the (77, k) plane resulting from a cut with the plane A^2 = 1 is also depicted. The plot 
evidences a region where A52 stays below 1. For comparison, in Figs. © and |2J), the corresponding 
dependence is shown for a state which only differs in the initial relative phase, 9 = 0° and 180°, 
respectively. The region A52 < 1 is not entered for such states. 

By keeping the effective time fixed at the value 77 = 77*, one can focus on some other dependencies 
as illustrated in Figs. Q and ©• In particular, Fig. (@J shows the dependence upon 9 and k. In 
this case, a reference contour plot is chosen at constant A52 = 0.8. The two contractivity islands 
characterized by (|33j) are clearly visible, and related to each other as in (|34j) . In Fig. ©, a similar 
plot depicts the dependences upon 5 and k. 
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B. Generalizations: Three-component cat states 

The above analysis shows that contractive quantum states, in the original Yuen's spirit, can 
be engineered through an appropriate choice of parameters in the class of cat states described 
by (|19|) , Several variants might be worth exploring in principle. For instance, cat states arising 
from the superposition of two Gaussians differing in their values of the parameters xq or Ao, or 
possessing non-zero initial momenta could be examined, as well as superpositions of non-Gaussian 
wave-packets. 

Another direction for generalizations is to consider multi-component cat states, for which con- 
tractivity could be enhanced via coherent interference effects between different pairs. We illustrate 
this possibility by focusing on the simplest generalization of the class (|19|) . leading to a family of 
cat states with three macroscopically distinguishable components. Thus, let us consider a wave- 
function of the form 

i^ssix) = N 3 {k + exp[-(x - x ) 2 /2A 2 ] + k exp[-x 2 /2A 2 } + k. exp[-(x + x ) 2 /2A 2 }} , (36) 

where the normalization constant is now given by 

|iV 3 | 2 = (7rAg)-V2 [\k + f + \k \ 2 + \k^\ 2 + (k* + k_ + k + kt)e- s2 

+ (k* + k + k + kl)e- 52 l A + {k*_k + fc-fcgje-**/ 4 ] _1 . (37) 

Here, S = xo/Aq as before, and we also continue to assume that the three Gaussians have zero 
initial momentum. Because a global phase factor is irrelevant, the state (|36|) is parametrized by 
five real parameters describing the relative amplitudes and phases, in addition to xq and Ao, which 
are taken as before to be positive. Similar to the two-component cat case, we arbitrarily choose one 
of the k coefficients to be real, thus setting ko = Co, k + = c + e td+ , fe_ = c~e %e ~ , with cq, c±,6± G R. 
With these conventions, the results for two-component cat states of the previous section can be 
recovered by letting Co = 0, 6 + = 0, and 6L = —8. We also define 

K = c\ + c\ + <?_ + 2c (c+ cos 9 + + c_ cos 0_)e" <52/4 + 2c+c_ cos(6>+ - 6L)e _<52 . (38) 

The expectations of the position and momentum observables on the state ()36|) take a more com- 
plicated form than in (|2T|) - (|22|) : 



(x)s3 = — c\ — c 2 _ + co(c+ cos 6+ — c- cos 6-)e s ^ xo , (39) 
K L 



/S3 



h co(c + sinl 



sin6»_)e-' 52 / 4 + 2c + c_ 



sin 



-5 2 



A 2 



c 2 _ + cq(c + cos + — c_ cos 8-)e 62 / 4 



\X)S3 



(40) 



Similarly, one also finds 



{xp)s3 = ^ttt? c + + c o + c - + 2co(c+ cos0+ + c_ cos^_)e 
2K 



-s 2 



+ 2c + c_cos(6» + -d-)e~ &2 - 2ic <5 2 (c+ sin6»+ + c_ sin6l_)e^ 2/4 



(41) 



Thus, unlike in the ordinary cat case leading to Eq. (|23|) . (xp) S3 acquires in general a non- vanishing 
real component. While one can verify that Im((xp) S3) = ih/2 and hence that a real value of the 
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correlation coefficient is correctly enforced, the full expression for A^ p 53(0) is less transparent than 
in the two-component case. In formal analogy to Eq. (|12j) . we may write 



A* S3 (0) = 2[Re({xp) S3 ) - (x) S3 (p} S 3 ) = -2K , (42) 



where 



C R 2 



- c (c+ sm9 + + c_ sm9_)Ke- s2/4 + (c (c+ sin0 + - c_ sm9-) e - s2/4 

+ 2c+c_ sin(#+ - 6-)e- 52 ^) (c\ - c 2 _ + c (c+ cos 9+ - c_ cos 0_)e _<52/4 ) . (43) 



The above function satisfies the property that 

C(c+ , c_ , 9 + , 9_ ) = C(c_ , c+ , ,0+ ) , (44) 

i.e., it is invariant under the exchange k+ *— > Thus, one may expect £ to be extremal for 
A; + = or c + = c_ and 9 + = 9-. The fact that under these conditions contractive behavior is 
possible for a cat state of the form (|3l)|) is easily verified by noting that (x)s3 = 0, (p)s3 = 0, and 
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C[ c+ =c_;0 + =e_] = -2 c c+ sin fl+e^ 2 / 4 . (45) 

Thus, ( > for 9 + G (180°, 360°), which in turn corresponds to a family of cat states with a 
negative correlation term. Note that this is in contrast with the two-component case, where a 
non-vanishing expectation of x (hence p) was found to be necessary for contractivity. 

The complete calculation of the position variance Ax| 3 (?7) as a function of the rescaled time 
is slightly simpler in the Schrodinger representation in this case. The essential steps and final 
equations are quoted in the appendix. Working as before in terms of a relative variance, 

As 3 (r?)= fi 3ir !\ =-+Bs + C 3 r l , (46) 

explicit expressions for the coefficients A3, 83,03 may be derived from the appendix. The be- 
havior of A53 has been analyzed numerically in the parameter space corresponding to the set 
i], k_, 9 + , 9-, 5, with = c+/c_ and k_ = co/c_. In this case, the global optimization leads 
to a minimum value 

A™"(r?*, <, 0*+, 01, S*) ^ 0.569 , (47) 

for optimal parameter values 77* ~ 1.270, k* + ~ 1.00, n*_ ~ 2.38, 9* + ~ 249°, 9*_ ~ 0* + , and 6* ~ 1.21. 
The fact that ~ 1.00 and 6*_ ~ 0+ within the numerical accuracy indicates that the minimum 
of A53 is indeed attained in the parameter regime expected from ([45 [I. A pictorial comparison 
between the behavior of the optimized contractivity figures of merit A^ 2 and A^ 3 as a function of r/ 
for two- and three-component cat states is displayed in figure ((SJ. Is it worth noting that not only 
A* S3 < A* S2 for T) > 0.5, but also the region where values lower than one are achieved appreciably 
widens for three-component contractive cat states. 

For an additional comparison between contractive and non-contractive behavior, it is also in- 
structive to directly examine the time dependence of the absolute variance for representative two-, 
three-component cat states, and for a Gaussian wave-packet with the same width parameter Ao- 
The results are shown in figure© for maximally contractive cat states. The monotonic increase of 
the variance for the Gaussian case, 

A4W = f (1 + 1 2 ). (48) 
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is clearly visible. The SQL behavior of Eq. Ijl]). which translates into A^Q L (r]) = A^tj, is reached 
at time rj = 1 in the non-contractive Gaussian case of (|48[) . The contractivity intervals for two- and 
three-component cat states are determined by the appropriate intersections of the variance curves 
A52, A53 with the SQL-line. Viewed in this way, the advantage of contractive states also manifests 
itself in the ability to preserve variances comparable to the initial Gaussian value at later times. 
In addition, the stronger contractivity of a three-component cat state is evidenced by the higher 
initial slope, implying both a deeper excursion below the SQL value and a longer contractivity 
interval. 

To summarize, the larger set of parameters available in more complex classes of Schrodinger 
cat states is capable of supporting stronger and better contractive properties. While further in- 
vestigation is needed to precisely characterize all the possible configurations, the fact that states 
with increased complexity may allow for a richer "self-interference" pattern (and hereby enhanced 
contractivity) is likely to occur beyond the specific case examined here. In particular, this is 
also somewhat reminiscent of the recent findings in 36], where cat states such as "compass- like" 
superpositions of four Gaussians are able to probe sub-Planck structures in the phase space. 



IV. IMPLICATIONS OF CONTRACTIVE CAT STATES 

The existence of contractive cat states immediately evokes, as in Yuen's original proposal, the 
potential for breaching the SQL of a free mass. However, a full assessment of such a potential 
requires a proper formulation of both the SQL itself and of various measurement-theoretical con- 
cepts such as quantum precision, resolution, and uncertainty. While addressing these issues goes 
beyond the scope of the present work, as far as the implications for the SQL are concerned we limit 
ourselves to a few preparatory remarks. 

As emphasized by Ozawa [2J|, because the SQL applies to repeated quantum measurements, 
one essential ingredient is the ability to perform a measurement using a single measuring appara- 
tus, the reading of which provides accurate information about the position of the free mass and 
simultaneously prepares it for the next position measurement, leaving it in a contractive state. 
A class of measurement schemes well suited for this task is offered by so-called Gordon-Louisell 
(GL) position measurementsM. The latter are naturally described within the general formalism 
of quantum operations 0, l38j| . Let a £ K be a real parameter, to be interpreted as the possible 
outcome of a position readout for a one-dimensional system, and let {|a}} denote the complete set 
of position eigenstates. If {^ a } is a family of normalized wave functions indexed by a, and p(0~) 
is the state of the system prior to the measurement, the state change and measurement statistics 
corresponding to a GL position measurement {| 1 I / a }(a|} are described by the following operation 
measure and effect measure, respectively: 

p(0+) = J daA a p(Q-)Al, A, = |*«)<o|,4 = |o>(*a|, (49) 



Prob(a G l\p(0~)) = Traced / daA\A a p(Q') 



Trace 



da \a) (a\p(0 



(50) 



for all Borel sets /CK, The distinctive feature of GL position measurement is that, according to 
(fi9")) . the posterior state p(0 + ) of the system is determined by \& a regardless of the prior state p(0~), 
the latter determining the outcome probabilities according to (|50j) . Because GL measurements are 
described by completely-positive quantum maps, every GL measurement is physically realizable in 
principle |24|. 
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GL contractive position measurements were originally invoked by Yuen (21)]. Given contractive 
states of the form @, the relevant family {^ajy may be constructed by fixing the parameters 
fjL,u,u, and by identifying xq = (x)y = a, for all real a. If we now look, for instance, at two- 
component cat states as in lfI5|). then a contractive state satisfying that (x)s2 = a can be found 
except for the (zero-measure) set I = {a = 0}. Let (3 > 1 a positive number. Then one may define 
a family {^ a }s2 a s follows: 



W I - ! ^> = ~ + 1) ; = 270°; xq = /fa, A = pa/8)} a>0, 

1 a/52 I {|«< = V(P + WW - 1); & = 90°; x = -pa, A = -Pa/5}} a<0. 1 j 

For fixed |a|, the two corresponding states in (|51|) have the same degree of contractivity thanks to 
the relationship ()34|) . By choosing 5 = 5* and adjusting /? in such a way that 1/k*, 
contractivity can then be brought close to optimality. Establishing whether a family of contractive 
cat states such as (|51|) actually supports a well-defined GL position measurement, and whether 
the latter may succeed in breakin g th e SQL, are interesting problems per se, which would call 
for an analysis along the lines of [24|. Once a successful strategy is identified in principle, an 
additional issue would be looking for an explicit interaction Hamiltonians implementing the desired 
contractive measurement scheme. For Yuen states, such a constructive problem was solved by 
Ozawa It is interesting to note that the proposed measurement model explicitly requires the 
initialization of the meter in a contractive state, and contractivity is subsequently transferred to 
the mass via a controlled interaction. In this context, cat states might then turn advantageous 
if preparation procedures simpler than the ones generating Yuen states could be exploited to 
appropriately initialize the meter. 

Independently of the SQL motivation, the possibility of experimentally demonstrating and char- 
acterizing contractive behavior for a free massive particle is both interesting as a fundamental 
quantum property, and for its possible implications in high-precision single quantum measure- 
ments. While being certainly challenging with present capabilities, accomplishing this goal could 
largely benefit from the extensive efforts which are under way to realize Schrodinger cat states in 
diverse physical systems. Starting from the pioneering proposals by Yurke and Stoler jjjjj 0] for 
the generation of optical cat states, and related following proposals 4l], 4^|, such efforts have led 



to the successful creation and manipulation of macroscopically distinguishable photon states in the 
cavity QED settin g [4j| , and have culminated in the generation of a mesoscopic Schrodinger cat for 
a trapped Be ion [4j|. Also, the creation of macroscopically distinguishable states has been re- 
cently reported for superconducting circuits 3 4^, 4fj], and proposals for entanglement of trapped 



electrons have been recently formulated 47l |. Interesting proposals exist for generating Schrodinger 
cat states of mechanical systems via controlled entanglement with microscopic degrees of freedom, 
including electron pairs in Cooper boxes j4^|, single photons in beam-splitter configurations |4fi| . 
as well as radiation-pressure controlled mirrors Finally, recent advances in 

ultracold atom physics open new perspectives for the possibility of crea ting Schrodinger cat states 
with ultracold atomic ensembles or gaseous Bose-Einstein condensates |5fiL 

Devices based on trapped ions or cold atoms seem especially promising for demonstrating and 
exploiting contractivity properties. In this case, if one of the strategies mentioned above actually 
succeeds for engineering a contractive Schrodinger cat state, trapping potentials could be switched 
off and the mass (ion or atomic cloud) be released to free evolution for a given time interval. The 
position variance could then be inferred from the spatial pattern of the fluorescent light emitted 
under the action of a light probe beam, and by repeating the experiment for different times of flight 
the relevant time dependence could be reconstructed. Besides demonstration, contractivity might 
find useful applications in high-precision atomic physics experiments. For instance, one possibility 
worth exploring in principle is trying to improve the accuracy of atomic-fountain clocks through 
proper refocusing of the atomic cloud on the fluorescence detection this could translate 
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into a decrease of the absolute uncertainty on the atomic populations [58| . A general observation 
which may be relevant in this context is that, because the behavior of the position variance would 
be unchanged for a mass subject to a linear potential, contractive states for a free mass would 
still exhibit contractivity under a uniform gravitational Geld (see also 0] for a related discussion). 
Thus, contractivity could in principle be exploited both in the vertical and horizontal directions. 

In spite of these promising avenues, a potentially serious concern arises from the fact that 
Schrodinger cat states may be especially fragile against the decoherence effects caused by the cou- 
pling to their surrounding environment , a feature which is regarded as crucial to understanding 
the quantum-to-classical transition 0, 0]. Under the assumption that the environment may be 
schematized as a harmonic bosonic bath, exact analyses are possible by adapting the results avail- 
able for the quantum Brownian motion mod el |ft l6(l| . or by directly resorting to results already 
obtained for a free damped quantum particle 61, 62]. However, it is also worth stressing that the 
decoherence behavior may be quite sensitive to various details of the system-environment coupling 
as well as environment properties, and there is hope that on one hand estimates based on general 
models might be in some cases overly pessimistic |fi3i ). and on the other hand decoherence effects 
might be effectively counteracted. In particular, stabilization procedure for cat states could be in 
principle designed based on the basis of both active-control schemes for the system and 
symmetrization schemes for the environment 66]. 



V. CONCLUSIONS 

We have shown that Schrodinger cat states provide a new class of states able to support con- 
tractive features in the sense originally proposed by Yuen in j^. Yuen contractive states lead 
in principle to an unlimited gain with respect to the SQL for a free particle but, apart from a 
proposal in the atomic physics context ji^], it is not clear how to prepare such states in practice. 
Contractive Schrodinger cat states, on the other hand, may enable one to capitalize on the intense 
effort which is ongoing to generate Schrodinger cat states with various quantum devices, ranging 
from trapped ions and atoms to superconducting circuits. 

Viewed from a broader perspective, our results enforce the conclusions reached in [jj, where the 
use of non-classical states and quantum resources of relevance for quantum information processing 
is also anticipated to improve our capabilities to gather information in high precision measure- 
ments. For Schrodinger cat states in particular, the identification of contractivity features provides 
independent additional support to the idea that highly delocalized states may develop enhanced 
quantum sensitivity as recently suggested by Zurek [36(. While several questions remain to be 
addressed in more depth, the present analysis thus points to a novel significance and potentially 
useful applications of Schrodinger cat states in the context of quantum- limited measurements. 



Acknowledgments 

L.V. gratefully acknowledges current support from the Los Alamos Office of the Director through 
a J. R. Oppenheimer Fellowship, and also wishes to thank the University of Padova for hospitality at 
the time in which these ideas originated. This work was also supported in part by Cofinanziamento 
MIUR protocollo MM02263577-001. 



13 



APPENDIX A: SCHRODINGER REPRESENTATION 



For completeness, we report the evaluation of the position variance in the Schrodinger repre- 
sentation. For two-component Schrodinger cat states, the starting point is the wave function 

ip S2 (x, 0) =N 2 {k + exp[-(x - x ) 2 /2A 2 ] + fc_ exp[-(x + x ) 2 /2A 2 ]} , 

where the normalization constant (also given in Eq. (|20|) ') is 



\N 2 



\k + \ 2 + |/t_| 2 + (k + k*_ + k*_k + )e 



-a 2 



and, as before, 5 = xo/Aq. The time evolution is obtained by solving the Schrodinger equation 
with effective time 77, 

dtp _ A 2 d 2 j> 
dr] 2 dx 2 ' 

which is easily accomplished by Fourier expanding the Gaussian states in terms of plane waves (see 
also |H3|). The resulting wave function may be written as: 



exp (— |atanry) 
(l + r? 2 ) 1 ^ 



(1 - irj)(x - x ) 2 \ 
k + exp ( „ „ ,„ — - s — I + 



k- exp 



2A (l + r ? 2 ) 
(1 - iij)(x + x ) 2 



2A (l + r ? 2 ) 

Thus, the required first and second momenta of the position observable are calculated as: 



and 



(^W^) = \N 2 \ 2 {irAl)^ 2 x 



(x\ c(x ^ = ^(ttA 2 ) 1 / 2 



(\k+\ 2 + \k- 



ir/(k + k*i — A;lfe+)e" 



: i + ? ? 2 ) A o , * 



-a 2 



,1 1 1 z ; / , I \^ 7 ? 2 )^0 2 2 \ -<5 2 

+ {k + k_ + k + k^) [ - rj x ) e 



The position variance Ax| 2 ( ? ?) 



)f S2 {x,ri) ' 



, \ agrees with the result obtained from 



the Heisenberg representation. 

In the case of three-component Schrodinger cat states as considered in Sect. IIIB, the initial 
wave function is 

ipssix, 0) = N 3 {k + exp[-(s - x Q ) 2 /2A 2 ] + k Q exp[-x 2 /2A 2 ] + k. exp[-(x + x ) 2 /2A 2 }} , 
with N3 (also quoted in Eq. I|37j)) given by 

|iV 3 | 2 = (vrAg)- 1 / 2 \\k + \ 2 + \k \ 2 + \k_\ 2 + (k* + k_ + k + k*_)e~ 52 
+ (k* + k + k + k*)e' s2 / 4 + (k*_k + k^e- 62 ^ 1 . 
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The time-evolved wave function ipsz (x, rj) can again be evaluated by solving the above Schrodinger 
equation. After cumbersome but straightforward calculations we obtain the average values of the 
position and its square as follows: 



<% S3 (^) = |iV 3 | 2 (7rA2)V2 Xo 



2 2 
c 1 — c_ 



+ ( kk*+k - k + kl)e- &2/A + i{k\k^ - k+k*_)e- 62 - Uk*_k Q - k^e^ 2 ^ 



and 



M + c o + ^)A'o(l + V*) + (4 + 



+ \{k%k + k + k* )e- 52 / 4 A 2 (l + ri 2 ) + \(k\k + k + kl)e- 52 ' A xl{l - rf) 



+ -(k* + k - k + k* )e- d ^ 4 x 2 r] + ^(k* + k_ + k+k*_ )e~ 52 A§(1 + ^)+ 
- (fc;A;_ + k+k*_ )e- s2 x 2 rj 2 + ^(k*_k + k_k* )e- s ' ''' A A 2 (l + V 2 ) 



+ -(k*_k + k^k* )e- &2 l A xl{l - if ) + Uk*_k - k_k* )e- 52 / 4 x 2 oV 



The position variance may be calculated from Ax 2 S3 (t]) = (^ 2 )^ S3 (x,r)) — (%) 2 p S3 (x ^ collecting 
terms of the same order in 77 and taking the ratio to the SQL value, the relative variance A^t/) 
may be cast in the form (|46|) , The above procedure was implemented numerically to calculate the 
optimal value A^ 3 and infer the corresponding parameters. 
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FIG. 1: The dependence of the contractivity figure of merit Ag2 upon the effective time 77 and k for optimal 
values 8* ~ 127° and 6* ~ 0.49. The intersection with the plane at constant A52 = 1 is also evidenced, 
showing the existence of an island where the variance of a Schrodinger cat state maintains values smaller 
than the SQL. 
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FIG. 2: The dependence of the contractivity figure of merit A$2 upon the effective time r\ and k for 9 = 0° 
and 5* ~ 0.49 as in Fig. 0J. No region with A52 < 1 is present in this case, therefore excluding any 
contractivity for the corresponding Schrodinger cat state. 




FIG. 3: As figure 0, but for 6 = 180°. 
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FIG. 4: The dependence of the contractivity figure of merit A52 upon 9 and k for optimal values 77* ~ 1.105 
and S* ~ 0.49. The two islands result from the section at constant Ag2 = 0.8, and are connected via the 
relationship given in equation 
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FIG. 6: Comparison of the optimized contractivity figures of merit A* S2 (solid blue) and A^ 3 (solid red) 
versus rescaled time for two- and three-component cat states as in equations and (j36|) . respectively. 
Except for the time, all the parameters are set to their optimal values as found upon minimization of the 
appropriate relative variance function. Note that lower values arc attained by A* S3 over a wider time interval 
than A* S2 . For reference, the unity value of a non-contractive state (dashed green) is also depicted. 



24 




Q L 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

FIG. 7: The dependence of the position variance Aa;| 2 (solid blue) and A| 3 (solid red) versus rescaled time 
for two- and three-component cat states. The variance is evaluated in units of Ag/2, which is the width 
of the probability distribution for each single Gaussian component of the cat. Except for the time, all the 
remaining parameters are set to the values maximizing contractivity as given in equations Q35f) and l|47l) . 
respectively. For comparison, the variance of a single-Gaussian state evolving as in equation (|48() from an 
initial variance Air 2 (0) = Aq/2 = 1/2 (dashed green) is also depicted in the same units. The SQL behavior 
corresponds, in these units, to a line (dash-dotted black) with slope 2 which is tangent to the Gaussian curve 
at r\ = 1 and intersects the contractive curves as shown. 



